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Abstract 

The isospin symmetry breakings of qq are investigated in the QCD sum 
rule method. The electromagnetic effects are evaluated following the procedure 
requiring that the electromagnetic effects for charged meson be gauge invari- 
ant. We find that the electromagnetic effects are also dominant in the isospin 
violations of the p mesons, which have been shown to be the case in the mass 
splittings of pions. The numerical results for the differences of pion decay con- 
stants and the masses of p mesons are presented, which are consistent with the 
data. 

PACS number(s): 13.40.Dk, 11.30.Hv, 12.38.Lg 
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The origin of mass differences in isospin multiplets has long been of great interest 
in nuclear and particle physics as a source of information about symmetry violations. 
Hadronic isospin violations are particularly important in that they arise from quark 
mass differences in Quantum Chromodynamics (QCD) (see Ref. fl[] for a review of 
the early work in this area), and of course electromagnetic effects. Among the first 
applications of the method of QCD sum rules was the study of isospin violations in 
the p — uj system^], in which it was recognized that the isospin splitting of the light- 
quark condensates can produce effects as large as the current-quark mass splittings 
and electromagnetic effects. Recently the QCD sum rule method has been used to 
study the neutron-proton mass difference p, f4j, the octet baryon mass splittings^ 
and the mass differences in the charmed meson systems (the D and D* scalar and 
vector mesons) 0. 

The isospin symmetry breakings of pions have been discussed extensively in the 
chiral perturbation theory, in which the electromagnetic effects are shown to be dom- 
inant in the mass splittings of pions. Thus, a consistent treatment of the electromag- 
netic effect is particularly important for the isospin symmetry breakings in the qq, 
q = u or d, systems. In comparison with hadronic quark models, the QCD sum rule 
method for calculating isospin splittings has the advantage that one can directly use 
QED field theory rather than rely on models to estimate Coulomb corrections. In 
Ref. P| , a systematic approach for the electromagnetic effects in the operator product 
expansion for the qq systems has been developed, in which the problem of the gauge 
invariance for the charged qq states was resolved. This enables us to investigate the 
violations of the isospin symmetry due to all three effects, the nonperturbative effects, 
the quark mass differences and the electromagnetic effects, consistently in the QCD 
sum rule method. The nonterturbative effect is due to the difference of the quark 
condensate, which arises from the up and down quark mass difference. The results 
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of our studies || in the heavy-light quark systems are in good agreement with the 
available data. The focus of this paper is to extend our early studies to the light-light 
quark systems, such as the 7r and p systems. 

The basic approach of the QCD sum rule is to study the two-point correlation 
function in the Wilson operator product expansion (OPE), defined by 

IV(g 2 ) = i J d A xe iqx < T{J^x) J v (0)) > 

= I% 2 )(^ - q 2 g flu ) + IliO 2 )^ , (1) 

where 

J»( x ) = ^( x ) l^{x) . (2) 
The IT(g 2 ) in Eq. can be written as 

Il( q 2 ) = U ( q 2 ) + U em (q 2 ) , (3) 

where ITo(g 2 ) is the leading term, and Ilem represents the contributions from the 
electromagnetic effects. The quantity n (g 2 ) has been evaluated in the literature. 
Thus we extend the results of n (g 2 ) in Refs. |2|, by including the flavor depen- 
dences of each term. In order to study the violations of the isospin symmetry, the 
electromagnetic effects n e m(<7 2 ) should also be written in the framework of the op- 
erator product expansion. The leading electromagnetic effects in this approach are 
the two-point functions from a two-loop perturbative contribution, whose Feynman 
diagrams are shown in Fig. 1. For the charge neutral current U.Q m (q 2 ), one could 
simply obtain the two-point functions by changing the gluons in QCD to the photons 
in QED in Fig. 1, since the two-point functions has been calculated in QCD|fJ, 0. 

n? m ( ? 2 ) = - i i 5 e^l„^, (4) 

where e q is the quark charge and \x is the infrared cutoff. Obviously, this result is 
gauge invariant. 
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Following the procedure in Ref . || , there are additional Feynman diagrams shown 
in Fig. 2 for the charged current so that the Ilg m (g 2 ) for the charged qq systems is 
gauge invariant. Thus, the two point function Hemio 2 ) for the charged qq systems is 

n em(<? 2 ) = -T^e gi e 32 -ln-|-. (5) 

The leading nonperturbative terms due to the electromagnetic effects are the four 
quark condensates, of which the Feynman diagrams are shown in Fig. 3 and Fig. 4. 
We present the expressions for IT(g 2 ) 

n p ± = -— (1 + — )ln— - — e u e d -\n — 

1 la 

+-(m u (dd) + m d (uu)) + — — (— G 2 ) 

12g 4 tv 

8 8 
+ 0^9s(uu)(dd) - —g 2 s ((uu) 2 + (dd) 2 ) 

+^e 2 e u e d ((uu) 2 + (dd) 2 ) + -L e 2 e 2 T (uu) (dd) , (6) 



and 



n o o = * (1 + ^) In - -4it^2 In 



P 47T 2 7T /i 2 47T 2 2 TV [I 2 

+^(m d (dd) + m u (uu)) + -^(— G 2 ) 



for the p systems, and 



n^i = —7-2(1 H ) In -— - —^e u e d - In — 

47T Z 7T yLT 47T^ 7T /i 
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+ -^(m u (dd) + m d (uu)) + -J-j^G 2 ) 
-a 2 s (uu)(dd) - -^- e 2 e u e d {(uu) 2 + (dd) 2 ), (8) 



and 



81g b ^ " ' 27q« 



47T 2 V 7T ^ 2 47T 2 2 7T /i 2 



1 la 

+—(m d (dd) + m u (uu)) + G 2 ) 

-^9s((^) 2 + (dd) 2 ) - ^e 2 (e 2 u (uu) 2 + e 2 d (dd) 2 ) (9) 

for the 7r systems, in which we have assumed the vacuum dominance and factoriza- 
tion hypothesis for the four quark condensates following Ref. 0. In the Appendix 
we present expressions in which the vacuum dominance approximation has not been 
implemented. The sensitivity to the possible deviation from vacuum dominance ap- 
proximation will be discussed later. 

The isospin symmetry breakings of these systems are determined by the difference 
between Il ± (g 2 ) — n°(g 2 ), and we have 



n p o - n p ± = - — -(e u - e d ) 2 -ln— %- + — — na(e u - e d ) 2 (uu) 2 
07r 7r pr 27q° 

4 „, . 1 



+ 7 ^g 2 ((uu) - (dd)) 2 + -{m u - m d )((uu) - (dd)) (10) 
9q b g 4 



for p states, and 



1 . x9 a, -cr 44 

' ) 

44 , , T _, 1 



n^o - n^i = - — -(e u - e d ) 2 - In —5- - — — -Txa(e u - e d ) 2 (uuf 
8ir z rc p z 27q° 



9t((uu) - (dd)) 2 + —(m u - m d )((uu) - (dd)) (11) 



for pions. The first two terms in Eqs. [10] and |IT] represents the electromagnetic effects, 
and the last two terms correspond to the isospin violations in the quark masses and 
condensates. Qualitatively, the contributions from the electromagnetic effects and the 
differences of the quark masses and condensates to the isospin symmetry breakings of 



7r and p states should have the same order of magnitudes. Thus, Eqs. [10] and [TT] show 
explicitly that the contributions from the differences of quark masses and condensates 
are of higher order relative to the electomagnetic effects. This has been shown in Ref. 



lfl for pions sometime ago, while Eq. 10 suggests that the electromagnetic effects is 



also dominant in the isospin symmetry breakings of p mesons. 



We adopt the standard one resonance plus a continuum model for the spectral 
density Imil p o(s) at the hadronic level: 



where f p o is related to the electronic width of the p meson and s p o pa 1.5GeV 2 is the 
p meson continuum threshold. We use ~ 0.75GeV 2 in the pseudo vector channel 
so that the sum rule is dominated by the pion and the ai(1260) meson contributes 
only to the continuum. Although it is impossible to extract the pion masses in the 
QCD sum rule framework, it is still possible to calculate their decay constants with 
the physical masses as the inputs [0, p]]. After making the Borel transformation and 
transferring the continuum contribution to the right hand side, the final sum rules 
are: 



Imn p o(s) = nf 2 5(s - m 2 ) + — (1 + + 



el + e 2 d a 

2 7T 



)0(s - S p o) 



(12) 




(13) 



and 




4tt 2 



1 



{(1 + 




Ml YlMl 
14 a s a 2 + a 2 _ 7 a ejaj + ejaj 
81 -k Mi 54 7T Mi 



(14) 



for p mesons, and 




4tt 2 



1 



{(1 + 




m u a d + m d a u (^G 2 ) 
Ml YlMl 
44a s a u a d 11 a a 2 + a 2 d 
81 7r M% 54 6u6 \ Mi ' 



(15) 
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and 
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B 



M| 12M| 
22 a s a 2 u + a 2 d 11 a e 2 u a 2 u + e 2 a 2 ^ 
+ 81 7T M| + 54tt M% 1 1 J 

where e T = e u - e d , e u = §, e d = -|, m u = 5.1MeV, m d = 8.9MeV, a 9 = -(27t) 2 (qq) , 

q=u, d, a = ? a s(-^i) = i n (iooAf a ) ' ^ s = anc ^ ^ ne va ^ ues °f the various 
condensates a g are standard@ a = = 0.55GeV 3 , and (^G 2 ) = 0.474GeV 4 . 

The mass sum rules for p mesons are obtained by taking the ratio of Eq. (|i~3 ) 



and the resulting equation of the first derivative with respect to ttt. The Am p 
is shown in Fig. 5. The threshold parameters for the stable sum rules are s p ± = 
1.502GeV 2 and s p o = 1.500GeV 2 . The mass of p ± is found to be 770 MeV, which 
reproduces the results in Refs. || [7j and is in good agreement with data. Our 
numerical analysis shows that contributions from the differences of quark masses and 
condensates are indeed negligible. We have the resulting Am p = —0.25 MeV, and 
this is not inconsistent with the experimental value from the Particle Data Group, of 
which Am p = 0.3 ± 2.2MeVQnj. 

With the physical pion masses m n ± = 139.57MeV and m n o = 134.98MeV as 
inputs, one could extract the decay constants of pions. The decay constant f n ± from 
these sum rules is 131 MeV, and we present the sum rule for Af n = f n ± — f n o in Fig. 
6. The corresponding threshold parameters are s n ± = 0.78GeV 2 and s T o = 0.70GeV 2 . 
Our result is Af n = 4.5MeV, while the data in the Particle Data Group for /„. suggest 
Af„ = 11.9 ±4.0 MeV|l3|. 



Our numerical analysis shows that the quantity Af n is dominated by the second 
term in Eq. [11], which comes from the four quark condensates. Thus, it is very 
sensitive to the vacuum dominance approximation adopted in our approach. The 



possible violation of this approximation was discussed in Ref. ||14j| . Following Ref. 

, we introduce a factor 1 + k in front of the QED induced four quark condensates 
based on the factorisation assumption in fll3|) , (|i~4l) , (|l~5|) and fllBl) . A nonvanishing k 
signals the breakdown of the factorisation hypothesis. We find that Am p is insensitive 
to k since it is dominated by the electromagnetic radiative correction. However, the 
A/tt is increased to 8.3MeV with k = 1, which is in good agreement with the PDG 
data. We estimate k = 0.5 ~ 1.5. More accurate data on Af n will determine how 
much the vacuum dominance approximation is being violated. 

In summary, we have presented a field theory calculation of the electromagnetic 
effects in the isospin symmetry breakings of the n and p mesons. We have shown 
explicitly that the isospin symmetry breaking induced by the elctromagnetic effects 
are dominant in the n and p systems. Our results show that the electromagnetic effects 
in isospin symmetry breakings of pions are much larger than that of p mesons. This 
is consistent with the available data, although we have not included the effects due to 
the p — ijj mixings in the mass splittings of p mesons. Certainly, more accurate data 
are needed to test our theoretical approach, and to determine the possible deviation 
from the factorization of the QED-induced four quark condensates. 

This work is supported in part by the Natural Science Foundation and the Doctoral 
Program of State Education Commission of China. The financial support by Peking 
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Appendix 



We present the expressions for the 7r° and p° systems in which the vacuum dominance 
approximation has not been implemented. 



All 2 7T H 2 47T 2 2 7T fl 2 

1 la 

-—(m d {dd) + m u (uu)) + 7^-7 (— G 2 



1 A a A a A a A a 

+^^ 2 [(^7M75yM)(«7' 1 75yM) + (d 7/ u7sy d)(d^ 5 —d)\ 

2 A a A a A a A a 

q=u,d,s ~ q=u,d,s 

+-^ e2 l e l(' a y^5u)( y uy^ 5 u) + e 2 d (d^^d)(d^^d)} 



2 

+ 7^e 2 [(e u M7 M M)( e q q^q) + (e d d^d)( ]T e g g 7 ^)], (17) 

q=u,d,s q=u,d,s 



for the p° system and 



47T 2 7r /i z 4n z 2 n /i z 

+\(m d (dd) + m u {uu)) + ~^{— G 2 ) 
q* IZq* 7T 

1 o r , A a w „A a , /T A a , wT „A a rtl 



2 A a A a A a A a 

+^<? 2 [(^y«)( E 97 M y9) + (rf7My rf )( E ^yg)] 
q=u,d,s ~ ~ q=u,d,s 

+ ^e 2 [e 2 (u 7 ^)(« 7 ^) + e 2 d (d^d)(d^d)} 

2 

+ ^-7e 2 [(e u M7^)( £ e^g) + (e d J 7 ^)( ]T e^g)] (18) 
q=u,d,s q=u,d,s 

for the 7T° system. Similarly we may derive expressions for il p ± and n^i. 
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Figure Captions 



1. The two loop perturbative corrections for charge neutral currents 

2. Additional diagrams generated by J^(x) for charged quark systems. 

3. The electromagnetic interaction induced four quark condensates for charge 
neutral currents. 

4. The additional four quark condensates for charged currents. 

5. The Borel mass dependence of Am p . 

6. The Borel mass dependence of Af n . 
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